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More than 30 years ago, Thouless introduced 
the concept of a topological charge pump [1] 
that would enable the robust transport of charge 
through an adiabatic cyclic evolution of the un¬ 
derlying Hamiltonian. In contrast to classical 
transport, the transported charge was shown to 
be quantized and purely determined by the topol¬ 
ogy of the pump cycle, making it robust to per¬ 
turbations [1, 2]. On a fundamental level, the 
quantized charge transport can be connected to 
a topological invariant, the Chern number, first 
introduced in the context of the integer quantum 
Hall effect [3, 4]. A Thouless quantum pump may 
therefore be regarded as a ’dynamical’ version of 
the integer quantum Hall effect. Here, we report 
on the realization of such a topological charge 
pump using ultracold bosonic atoms that form a 
Mott insulator in a dynamically controlled optical 
superlattice potential. By taking in-situ images of 
the atom cloud, we observe a quantized deflection 
per pump cycle. We reveal the genuine quantum 
nature of the pump by showing that, in contrast 
to ground state particles, a counterintuitive re¬ 
versed deflection occurs when particles are pre¬ 
pared in the first excited band. Furthermore, we 
were able to directly demonstrate that the system 
undergoes a controlled topological phase transi¬ 
tion in higher bands when tuning the superlattice 
parameters. 

Charge pumping in solid-state systems has received 
much attention, mainly due to its potential for realiz¬ 
ing novel current standards [5, 6], but also for charac¬ 
terizing many-body systems [7, 8]. Quantized transport 
of electrons without an external bias was observed in 
tunnel junctions with modulated gate voltages [9], in 
ID channels using surface acoustic waves [10] and in 
quantum dots [11]. While the latter was used to re¬ 
alize an adiabatic quantum pump [12, 13], topological 
charge pumping has so far remained out of reach in con¬ 
densed matter experiments. In engineered bosonic sys¬ 
tems, non-quantized topological pumping of edge states 
was observed in ID quasicrystalline photonic waveguide 
arrays [14, 15]. In cold atomic systems, topological Bloch 
bands have been realized, ranging from the Su-Schrieffer- 


Heeger (SSH) model [16], the Hofstadter model in real 
space [17, 18] and synthetic dimensions [19, 20] to the 
Haldane model [21] . Their geometric features have been 
probed using novel interferometric [16, 22] and transport 
probes [21, 23] that have e.g. enabled the direct mea¬ 
surement of the Chern number through bulk topological 
currents [23] . 

Due to their versatility, ultracold atoms in optical su¬ 
perlattices constitute an ideal system for the implemen¬ 
tation of quantized topological charge pumps [24-26] . 
A superlattice is formed by superimposing two lattices 
with periodicities di and ds = adi, a < 1. This creates 
a potential 14 sin^ {-Kx/dg 7r/2) -|- Vi sin^ {nx/di — 7>/2), 
where 14 (M) denotes the depth of the short (long) lat¬ 
tice, respectively. The relative position of the lattices 
is determined by the variable phase ip. When applying 
Bloch’s theorem, the resulting single-particle Hamilto¬ 
nian 'H{kx,7>) is periodic in both the quasi-momentum 
kx £ ]—7r/di, 7r/di] and the superlattice phase (p. 

A cyclic pumping scheme can be realized by adia- 
batically varying p where one cycle corresponds to a 
change of 27r, i.e. moving the long lattice by d\. A 
particle which is initially in an eigenstate \un{kx, po)) 
of 'H{kx,Po) in th® n-th band follows the correspond¬ 
ing instantaneous eigenstate \un{kx, p))■ Even with per¬ 
fect adiabaticity, however, the particle acquires a small 
admixture of states from other bands, proportional to 
dtp. This gives rise to an anomalous velocity in = 
rindtp, determined by the Berry curvature Fln{kx,p) = 
(^(^dipUn\dk„.ttn'} (,dj^^Un\d(p'iin'}') [27, 28[. Hence, chang¬ 
ing p induces a motion of the particle which, depending 
on the sign of is either in the same or opposite direc¬ 
tion as the moving lattice. The displacement after one 
cycle is obtained by integrating in and can in principle 
take any arbitrary value. 

For a filled or homogeneously populated band, how¬ 
ever, the average displacement of the entire cloud per 
cycle can be related to a 2D topological invariant, the 
Chern number of the pumping process: 

difdkx (1) 

27 ^ Jfbz Jo 

Here, FBZ denotes the first Brillouin zone of the super¬ 
lattice. Such a case can be realized with both fermions. 
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by placing the Fermi energy in a band gap, and bosons, 
by preparing a Mott insulator in the n-th band as in our 
experiment. During one cycle, the cloud’s center-of-mass 
(COM) position changes by Vnd-i and is quantized in units 
of d\ because Vn can only take integer values. Unlike in 
classical systems, this motion can be faster {vn > 1) or 
even opposite < 0) compared to the one of the under¬ 
lying lattice [29]. As the displacement is proportional to 
a topological invariant, it neither depends on the pump¬ 
ing speed, provided adiabaticity still holds, nor on the 
specific lattice parameters as long as band crossings do 
not occur. Hence, the transport is highly robust against 
perturbations such as interaction effects in fermionic sys¬ 
tems or disorder [2]. 

The connection to a Chern topological invariant shows 
that charge pumping in a ID superlattice is closely re¬ 
lated to the integer quantum Hall effect, where charged 
particles move in 2D in the presence of a perpendicular 
magnetic field. Indeed, the adiabatic variation of ip is 
equivalent to the threading of magnetic flux through a 
cylinder [1] that generates an electric field in the orthog¬ 
onal direction and leads to the quantized Hall conduc¬ 
tance [4]. A direct analogy between these systems can 
be made in two limiting cases: the quantum sliding lat¬ 
tice with Us —>■ 0 and the deep tight-binding limit with 
Us 3> /{AEt- s)- The former matches with the free par¬ 

ticle limit of Landau levels that are associated with Chern 
numbers Vn = -\-\. In this limit, one can understand the 
pumping by following Laughlin’s argument [30], where 
the threading of magnetic flux leads to a sliding of the 
localized oscillator centers. In the latter, one obtains 
the generalized ID Harper model [31, 32] which, using 
dimensional extension [33, 34], can be mapped onto the 
2D Harper-Hofstadter-Hatsugai (HHH) model describing 
non-interacting particles on a 2D square lattice with a 
uniform magnetic flux 27ra and nearest as well as next- 
nearest-neighbor hopping [31, 35-37] (see Supplementary 
Material). In this mapping, ip corresponds to the trans¬ 
verse quasi-momentum ky. In the following, unless men¬ 
tioned otherwise, we will focus on the Wannier tunnel¬ 
ing regime U ^ Vx /where - even though the 
direct mapping to the ID Harper model breaks down - 
the pump is characterized by the same Chern number 
distribution. 

In the experiment, we use a superlattice with di = 2ds 
and two sites per unit cell (Fig. la). In the tight-binding 
limit, it is described by the Rice-Mele Hamiltonian [38] 

H{p)=-J2 {Mv’)blnam + M^)aln+ibm + h-c.^ 

m 

m 

(2) 

where dh^ (dm) ^-nd (bm) the creation (annihila¬ 
tion) operators acting on the even (left) and odd (right) 
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Figure 1. Topological charge pumping in an optical super lat¬ 
tice. (a) Superlattice potential created by superimposing two 
lattices with periodicities di and da = di/2. Here, Ji and J 2 
denote the tunnel couplings and A the energy offset between 
neighboring sites, (b) Evolution of the ground state Wannier 
function (red) during the first half of a pumping cycle. The 
long lattice is shifted to the right when increasing the super¬ 
lattice phase p. An atom initially localized in a symmetric 
superposition on a double well tunnels to the lower lying site 
and thereby follows the motion of the long lattice in a quan¬ 
tized fashion. A classical particle, on the other hand, stays 
at a fixed position as the individual sites do not move, (c) 
Pumping cycle in the (Ji — J 2 )-A parameter space. Vary¬ 
ing p from 0 to 27r corresponds to a closed path around the 
degeneracy point at the origin where A = 0 and Ji = J 2 . 
(d) Berry curvature fli of the lowest band as a function of 
p and the quasi-momentum for a lattice with U, = 10i5r,a 
and Vi = 20i?r,i, where Er,i = /i^/(2maAf) denotes the cor¬ 
responding recoil energy, rria, the mass of an atom and \i the 
respective wavelength. The panel below shows the Berry cur¬ 
vature averaged along k^, as seen by a particle localized in a 
single double well, which is peaked at the symmetric double 
well configurations p = ln,l G 1.. The graph on the right 
shows Qi{ka:) for p = 0. 


site of the m-th unit cell, respectively. Ji, J 2 denote the 
tunnel couplings and A the energy offset between neigh¬ 
boring sites. For A = 0, this Hamiltonian reduces to the 
SSH model [39]. 

The mechanism underlying the pumping can also be 
understood on a microscopic level. Shifting the phase p 
changes the shape of the potential (Fig. lb) and modifies 
Ji, J 2 and A periodically (Fig. Ic). At :/? = 0 (A = 0, 
Ji > J 2 ), a ground-state particle localized in a single 
double well is initially in a symmetric superposition of 
residing on both the left and right sites. With increasing 
p, i.e. shifting the long lattice to the right, the double 
wells are tilted (A > 0) and the atom tunnels to the lower 
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Figure 2. Center-of-mass (COM) position of the atom cloud as a function of the pumping parameter ip for the lowest band 
with 14 = 10.0(3)-Er.B and Vi = 20(1) (a) Detailed evolution of the measured COM position during a pump cycle. The 

step-like motion is caused by tunneling of the atoms to the lower lying sites. The solid black line depicts the calculated COM 
motion of a localized Wannier function. Each point is the average of ten data sets and the error bar shows the error of the 
mean. COM positions are determined differentially comparing a sequence with pumping to a reference sequence of the same 
length, but with constant phase (p = 0. One data set is obtained by averaging ten images each, taken in alternating order, and 
subtracting the resulting COM positions. The inset illustrates the motion of a localized Wannier function in the first band 
during a full pump cycle, (b) COM displacement and site populations for multiple pump cycles in the positive and negative 
pumping direction. The COM positions are averaged over ten data sets and the error bars depict the standard deviation. The 
dashed black line shows the ideal motion of a localized Wannier function and the gray line is a fit of the Wannier COM to 
the data assuming a finite pumping efficiency of 97.9(2)% per each half of a pump cycle (see Supplementary Material). The 
inset displays the population imbalance between even and odd sites as a function of ip with rie (wo) the fraction of atoms on 
even (odd) sites. Each data point is the average of five measurements and the error bars indicate the corresponding standard 
deviation. The gray line is obtained by fitting the calculated even-odd distribution of the Wannier function using the same 
model as for the COM displacement which yields an efficiency of 98.7(1)% leading to a slight decrease of the imbalance over 
time. 


lying site on the right. The sign of Ji — J 2 is reversed 
at If = 7r/2, where the tilt is largest, and at ip = tt the 
lattice forms symmetric double wells again, but shifted 
by one short lattice constant to the right. The atom, 
which remained on the lower site for large A, delocalizes 
over the shifted double well as A becomes comparable 
to J 2 and therefore has moved by di/2 during the first 
half of the pumping cycle. In the second half, the same 
procedure is repeated, but shifted by one site. After one 
cycle, the lattice configuration is identical to the starting 
point, but the atom ends up in the double well next to 
the initial one. In contrast to this, a classical particle 
would not move because the positions of the individual 
sites do not change. This illustrates the importance of 
quantum tunneling for the pumping. 

During the pump cycle, the system moves along a 
closed trajectory in the (Ji — T 2 )-A parameter space. It 
encircles the degeneracy point at A = 0 and Ji = J 2 , 
where the two bands of the Rice-Mele model touch, 
and smoothly connects the topologically distinct phases 
Ji < J 2 and Ji > J 2 of the SSH model (Fig. Ic). The 
Berry curvature and thus the motion of the atoms is 
peaked around p = hr, I € Z where the tilt changes sign 
and the atoms tunnel to the neighboring sites (Fig. Id). 


The Chern numbers of the pump cycles in the two bands 
of this model are = +1 and 1 Z 2 = —1, giving the 
same Chern number distribution as the HHH model with 
a = 1/2 [37]. Their sum vanishes because these bands 
emerge from the topologically trivial lowest Bloch band 
of the short lattice. 

The experimental sequence starts by preparing an 
n = 1/2 Mott insulator of ®^Rb atoms in a 3D optical 
lattice with at most one atom per unit cell in the ground 
state of symmetric double wells with = 0 and Ji 3> J 2 
(see Methods). Due to the large on-site interaction, each 
atom is localized on an individual double well, resulting 
in a homogeneous delocalization over the entire first Bril- 
louin zone. The pumping is performed by adiabatically 
shifting the phase (f of the long lattice (see Methods) and 
the resulting motion of the atoms is tracked by measuring 
the COM position of the cloud in situ. The displacement 
during one cycle is indeed quantized and occurs in steps 
(Fig. 2a) - unlike the underlying linear motion of the 
long lattice. The cloud moves by one lattice constant d\ 
per cycle as expected for = +1 and the steps appear 
around ip> = /tt, / € Z, where the atoms tunnel from one 
side of the double wells to the other. When perform¬ 
ing multiple cycles, the cloud keeps moving to the right 
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Figure 3. Transition from a quantum sliding lattice to the 
Wannier tunneling limit for the lowest band. Differential 
deflection Aa; after one pump cycle for positive and nega¬ 
tive pumping direction for various lattice depths 14 at Vi = 
25(1) Er.i. Each point consists of ten data sets comparing 
the COM position of ten averaged images for both directions. 
The error bars depict the error of the mean. For the data 
points in the tight-binding regime, the insets show the cor¬ 
responding pump cycles in the (Ji — J 2 )-A parameter space. 
For 14 = SEr.s, the two-band model breaks down for large 
tilts such that Ji , J 2 and A are not well-defined. The dashed 
line therefore connects the points where the gap between the 
second and third band becomes smaller than lOJi for (p = 0. 


whereas it propagates in the opposite direction for the 
reversed pumping direction (/? < 0 (Fig. 2b). The small 
deviation from the expected displacement for the motion 
of ideal Wannier functions can be attributed to a finite 
pumping efficiency due to non-adiabatic band transitions 
and the additional trapping potential. 

The step-like transport behavior can also be observed 
in site-resolved band mapping measurements (inset of 
Fig. 2b) which determine the number of atoms on even 
and odd sites. As for the COM position, a step occurs 
in the even-odd distribution whenever a symmetric dou¬ 
ble well configuration is crossed at (p = Itt,1 G Z. Using 
the measured even-odd fractions, one can estimate the 
transfer efficiency, i.e. the fraction of atoms transferred 
from site * to i -b 1 at each step. This is equivalent to 
the fraction staying in the lowest band during one half of 
the pumping cycle and allows to quantify the adiabatic- 
ity of the pumping protocol. From our data we obtain 
an efficiency of 98.7(1)% (see Supplementary Material). 
With the same model, the ideal COM displacement can 
be fitted to the measured positions which yields an ef¬ 
ficiency of 97.9(2)% with the small additional reduction 
most likely being caused by the trap. 

Due to the topological nature of the pumping, the dis- 
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Figure 4. Cloud displacement and site occupations for the 
first excited band with 14 = 10.0(3) and Vi = 20(1) Er,i. 
The main plot shows the evolution of the COM position for 
up to 1.5 pump cycles. The data points are averaged over ten 
data sets with the errors being the error of the mean. The 
dashed black line indicates the motion of a localized Wannier 
function of the second band and the gray line is a fit to the 
data using the same model as in Fig. 2b, giving a transfer 
efficiency of 97(2)%. The inset illustrates the evolution of the 
Wannier functions of the first and second band during a pump 
cycle, showing a deflection in opposite directions. The lower 
plot shows the imbalance of the fraction of atoms on even (rie) 
and odd sites (rio) averaged over 3-6 measurements each with 
the error bars depicting the standard deviation. The gray 
line is a fit of the even-odd distribution of the corresponding 
Wannier functions with a pumping efficiency of 96.7(3)% (see 
Supplementary Material). 


placement per cycle for the lowest band does not depend 
on the path in the (Ji — J 2 )-A plane as long as it encom¬ 
passes the degeneracy point. Moreover, it is independent 
of 14 since the sliding lattice and the tight-binding Thou- 
less pump are topologically equivalent for the first band 
and connected by a smooth crossover without closing the 
gap to the second band. To verify this, we measured the 
deflection of the cloud with V\ = 25(1) Ej- i for various val¬ 
ues of Us- For all parameters, the resulting displacements 
are consistent within the error bars (Fig. 3). 

The excited band in the Rice-Mele model exhibits 
counter-propagating charge pumping with V 2 = —1, i.e. 
the atoms are expected to move in the opposite direction 
as the long lattice. This underlines the pump’s intrinsic 
quantum mechanical character as such a motion could 
not occur for classical particles. To study this, the atoms 
were prepared in the second band at tp = 0 (see Methods) 
and pumping was performed with identical parameters 
as for the first band in Fig. 2. When moving the lattice 
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to the right {ip > 0), the cloud indeed shifts to the left 
with X < 0 (Fig. 4). This is further confirmed by the 
measured site occupations, showing that the behavior is 
exactly reversed compared to the lowest band. While 
the first band is localized on the lower site for large A, 
atoms in the second band are found on the upper site and 
are therefore transported in the opposite direction. The 
slightly larger deviation from theory is mostly due to the 
finite lifetime of atoms in the higher band. This leads to 
a lower transfer efficiency of 96.7(3)% obtained from the 
site populations, in agreement with the value from the 
fit to the COM positions of 97(2)% (see Supplementary 
Material). 

By varying Vs/V^, one can study the transition be¬ 
tween the topologically very different regimes of the slid¬ 
ing long lattice with 14 = 0 and the Wannier tunnel¬ 
ing limit for 14 > /(AE-c^s)- For the latter, the Chern 
numbers of the bands alternate between V 2 n+i = +1 and 
V 2 n = — 1 for d\ = 2ds, like in the HHH model with a 
flux of TT. This causes the opposite deflections for the 
first and second band observed in Fig. 2 and Fig. 4. In 
the other limit 14 —>■ OE^-^s, however, each band corre¬ 
sponds to a single Landau level with i/n = -|-1. Between 
these two limiting cases, as 14 is decreased, an infinite 
series of topological phase transitions occurs where two 
bands touch and exchange their Chern numbers (Fig. 5a). 
Thereby the negative Vn are successively transferred to 
higher bands and they become positive for the lower 
bands. 

While the deflection in the first band is independent 
of 14, there is a transition for the second band where the 
gap to the third band closes and 1/2 changes from —1 to 
-|-1. This transition can be mapped out by measuring 
the displacement of the cloud as a function of the lattice 
depths 14 and Vi (Fig. 5b). The direction of the motion 
reverses when crossing the transition, which occurs at 
14 = /{AEj-^s) in the tight-binding limit. For small 14, 

atoms in the second band move in the same direction as 
the lattice, whereas they move in the opposite directionin 
the Wannier tunneling regime. 

In conclusion, we have demonstrated the implementa¬ 
tion of a topological charge pump using ultracold atoms. 
Unlike previous experiments studying non-quantized 
pumping of edge states, we observe a quantized response 
in the bulk. Combining such a pump with novel tech¬ 
niques for engineering optical potentials at the single-site 
level [40, 41] would allow for a direct observation of 
edge states in finite systems [42] and to study their 
transport properties [14] . Furthermore, by adjusting the 
ratio of the lattice constants a = ds/d\, one can realize 
a wide variety of commensurate and incommensurate 
superlattices where the Chern number of the lowest 
band can in principle take arbitrary integer values. In 
addition to the negative deflection for < 0 shown 
here, another counterintuitive case of charge pumping 
can occur in these systems where the atoms move faster 
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Figure 5. Topological phase transition in the first excited 
band, (a) Band structure aX ip = O.Stt, where the band cross¬ 
ings occur, and Chern number distribution versus 14 with 
V\ = 25i5r,i. The shaded areas illustrate the width of the 
Bloch bands, which were calculated by numerical diagonal- 
ization of 'H{kx,p)- The color indicates the Chern number 
of the corresponding pump cycle with red being -|-1 and blue 
— 1. Starting from alternating Chern numbers in the Wan¬ 
nier tunneling limit for large 14, the negative Chern numbers 
successively propagate towards higher bands in a series of 
topological phase transitions when lowering 14, giving rise to 
a uniform distribution with Vn = +1 in the Landau limit of a 
sliding long lattice. The circle highlights the crossing of the 
2nd and 3rd band leading to the topological transition studied 
in b. (b) Differential deflection Aa: between single pump cy¬ 
cles in opposite directions for the second band as a function of 
the lattice depths 14 and Vj- At small 14, V 2 = -fl and atoms 
move in the same direction as the lattice, but as 14 increases 
a topological transition occurs at 14 = Vj^/(4i5r,s), indicated 
by the black line, where 1/2 suddenly changes to —1 and the 
direction of motion is reversed. Each square corresponds to 
one data set averaged over 10-12 pairs of images. 


than the long lattice for > 1 [29] . By adding a spin 
degree of freedom, the pumping scheme can be used to 
implement the Z 2 spin pump [43, 44] in a spin-dependent 
super lattice [45] . Moreover, extending the Thouless 
pump to 2D systems would enable the realization of an 
analogon of the 4D integer quantum Hall effect [34, 46] . 


Note Recently, we became aware of similar work by 
S. Nakajima et al. [47] implementing the Thouless pump 
with fermionic atoms. 


We acknowledge insightful discussions with F. Grusdt 
and S. Kohler. This work was supported by NIM and the 
EU (UQUAM, SIQS). M. L. was additionally supported 
by ExQM and O. Z. by the Swiss National Science Foun¬ 
dation. 
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METHODS 


Initial state preparation in the first band All se¬ 
quences started by loading an n = 1 Mott insulator of about 
3000 ®^Rb atoms in the lowest band of a 3D optical lattice. 
The lattice was created by three orthogonal standing waves 
with wavelengths Ai = 1534 nm for the long lattice along the 
2;-direction and Xy = 767 nm and Az = 844 nm along the y- 
and 2-axes, respectively. They were ramped up in 150 ms to 
a depth of Vi = 30(1) Er,i,i € {l,y,z}. The superlattice po¬ 
tential was created by adding another lattice with As = Ai/2 
along the ^-direction. To prepare the atoms in the lowest 
band of the superlattice, the lattice sites along the a:-axis 
were split symmetrically with ip = 0.00(l)7r by ramping up 
the short lattice to 14 = 10.0(3) i?r,s within 10 ms while si¬ 
multaneously lowering VI to 20(1) such that Ji J2. It 
was verified with direct band mapping measurements that the 
atoms are homogeneously distributed over the entire band. 

Initial state preparation in the second band For the 

preparation in the second band, the splitting was instead per¬ 
formed at 95 = 0.11(l)7r and 14 was ramped to 30(1) Er,s- The 
phase was then changed non-adiabatically to ip = —0.08(l)7r 
in 20 ms to transfer all atoms to the excited band. After that, 
the lattices along x were lowered to 14 = 10.0(3) Er.s and 14 
to 20(1) respectively, in 2 ms and the phase was moved to 
ip = 0.00(l)7r in 10 ms. This brings the atoms into the excited 
state of symmetric double wells with almost perfect efficiency. 

Sequence for pumping The pumping cycle was imple¬ 
mented experimentally by slightly changing the laser fre¬ 
quency and thereby shifting the phase ip of the long lattice 
along the a:-direction. For this, two separate lasers were used 
with one laser covering the range from —0.50(l)7r to 0.62(l)7r 
and the second laser from 0.62(l)7r to 1.50(l)7r (see Supple¬ 
mentary Material). 
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MAPPING TO THE HOFSTADTER MODEL IN 
THE DEEP TIGHT-BINDING REGIME 

In the tight-binding approximation, the dynamics of 
ultracold atoms in an optical superlattice potential with 
arbitrary lattice constants di and ds = adi, a < 1 can be 
described by the following Hamiltonian: 

((-^0 + SJ^ip)) aL+iOm + h-c.) 

A / ^ ^mip)^7n^m 

m 

Here, d\^ {d^) is the creation (annihilation) operator act¬ 
ing on the m-th lattice site and J denotes the hopping 
amplitude between neighboring sites in the short lattice. 
5Jm is the modulation of the hopping amplitude and 
Am the on-site energy of the m-th site, both of which 
are induced by the long lattice and depend on p. For 
14 ^ s), Eq. S.l reduces to a generalized ID 

Harper model [SI, S2] with 



Figure SI. Square lattice in the presence of a uniform mag¬ 
netic field with nearest-neighbor hopping Jx = Jo, Jy = A/A 
and next-nearest-neighbor hopping Jd = 6J/A using the gauge 
of Eq. S.4. 


Expanding dl^ and ^ into their Fourier components, 
aln,v = T,n and gives 


6Jm{p) = cos(27rQ:m — p) 

A , , 

Am{p) = cos(27ra(TO - 1/2) - p) 

Following the approach of dimensional extension [S3], 
the periodic parameter p can be regarded as an extra 
dimension, i.e. a realization for a given p corresponds to 
a ID subset of a 2D system. The pumping then samples 
the full 2D parameter space which is the reason why the 
displacement can be related to a 2D topological invari¬ 
ant. The 2D Hamiltonian onto which the ID superlat¬ 
tice maps in this limit can be found by thinking of the 
ID Hamiltonian as a single Fourier component of the 2D 
one and performing an inverse Fourier transform. For 
this, we relabel the creation and annihilation operators 
as and dm and define a 2D Hamiltonian 


i72D = ^^ H{p)dp (S.3) 


H 2 T 0 = — ^ Jod. 


t z 

'm+l,n m,n 


A 




,+22-7ra(m—1/2) -1 


+227ram;^t 

A ^ 


dJ 


E OJ 

A 


— i27rci:m 


+ h.c. 


(S.4) 


which is precisely the 2D Harper-Hofstadter-Hatsugai 
model describing non-interacting particles on a 2D square 
lattice in the presence of a uniform magnetic flux 27rQ! per 
square plaquette with nearest and next-nearest neighbor 
tunneling as illustrated in Fig. SI [S4[. Under this map¬ 
ping, the superlattice phase p corresponds to the trans¬ 
verse quasi-momentum ky and in this sense the pumping 
protocol is equivalent to performing Bloch oscillations in 



the 2D lattice by applying a gradient. For our particular 
case of ds = di/2, one obtains a flux of tt per plaquette. 

MAPPING OF THE SLIDING LATTIGE TO 
LANDAU LEVELS 

In the opposite limit of the sliding lattice with I 4 —>■ 0, 
a similar analogy can be made to the Landau levels of 
a free particle in an external magnetic held in 2D. The 
Hamiltonian of a single particle in the long lattice reads 

Hi = ——h Vi sin^ {irax/ds — ^/2) (S.5) 

2ma 

For a deep lattice and when only considering the lowest 
bands, the lattice potential can be expanded in orders of 
X — Xm around each lattice site Xm (he. the minimum of 
the potential). In second order this gives: 



0 160 ^ 260 ^ 360 


t (ms) 

Figure S2. Schematic of the experimental sequence with the 
lattice and phase ramps for the initial state preparation as well 
as for two pumping cycles. The light and dark lines indicate 
the two different lasers used sequentially to overcome tuning 
range limitations. The phase ip is shown modulo 2n. 


y.O 2 2 

rr pi . 


(S.6) 


Neglecting the coupling between neighboring sites, this 
splits into a series of decoupled Hamiltonians for each 
lattice site, each of which is a ID harmonic oscillator. 

On the other hand, the motion of a particle with charge 
q in an external magnetic held B = Be^, along the z- 
direction is described by the Hamiltonian 


i^LL = 


2m, 

pI 




2ma 2m, 


{Py - qBxf 


(S.7) 


with the vector potential A = BxSy given in the Landau 
gauge. As the above Hamiltonian commutes with the 
operator for the transverse momentum py, one can define 
a common set of eigenstates of both Hll and py. For a 
given state with momentum hky, Eq. S.7 can also be 
rewritten as a ID harmonic oscillator 


TOPOLOGICAL PHASE TRANSITION 

For ds = d\/2, the phase transition to the Wannier tun¬ 
neling regime, where the pump is characterized by the 
same Chern number distribution as the 2D HHH model, 
occurs when the second and third band cross at the stag¬ 
gered configuration p = 7r/2 during the pumping cycle. 
In the tight binding limit, the crossing point can be deter¬ 
mined by comparing the energy of the first excited state 
on the lower site with the one of the ground state on 
the higher site. For a sufficiently deep short lattice, they 
can be approximated by harmonic oscillator states with 
En = {n -\- 1/2 )}vjj. The on-site trapping frequency uj is 
given by 


^ m^dl ~ 

At tp = 7r/2, the tilt is largest with A = and the bands 
touch if A = fiw. Hence, the transition occurs at 

Vi = 2,/VsWys (S.IO) 


Hll = 


2m, 




Hkv 


maWr 


(S.8) 


where Wc = qB/m^, is the cyclotron frequency. Compar¬ 
ing this with Eq. S.6, one can see that the state in the 
n-th band localized at the lattice site Xm corresponds to 
the state in the n-th Landau level with a transverse mo¬ 
mentum hky = m^ujcXm- As in the tight-binding regime, 
the number of magnetic flux quanta per unit cell of the 
(non-existent) short lattice is given by the ratio of the 
lattice constants a = ds/d\ in this mapping. 


A similar derivation can be made for ds = d\/[2j),j € 
N which gives the same transition point. For all other 
ratios, Eq. S.IO gives a lower bound of k)(I4) and the 
transition will in general take place at slightly larger Vy. 


EXPERIMENTAL SEQUENCE 

Figure S2 illustrates the detailed experimental se¬ 
quence for the initial state preparation of the atoms in 
the lowest band of the superlattice potential and two ex¬ 
emplary pump cycles. As discussed in the Methods of 
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the main text, the sequence started by loading an n = 1 
Mott insulator in 150 ms in a 3D optical lattice created 
by the long lattice and two orthogonal standing waves 
with depth Vi = 30(l)i?r,iD € Subsequently, 

the long lattice sites were split symmetrically by ramp¬ 
ing up the overlapped short lattice with As = Ai/2 to 
14 = 10.0(3) i?r,s at = 0- Simultaneously, the long 
lattice was lowered to Vi = 20(1) such that Ji J 2 . 

For the pumping, the phase of the long lattice was 
moved by changing the laser frequency slightly. Due to 
the limited tuning range of a single laser, a successive 
hand over between two independent laser systems (in¬ 
dicated with light and dark lines in Fig. S2 for the dif¬ 
ferent lasers) was implemented in order to extend the 
pumping range to arbitrary values. Each pump cycle 
consists of four segments with s-shaped phase ramps 
ip € [0, 0.627r], [0.627r, tt], [tt, I.Stt] and [1.57r,27r] of 50ms 
duration to minimize the probability for non-adiabatic 
transitions to higher bands at the symmetric configura¬ 
tions (p = hr, I G J.. For the measurements in Fig. 3 
and Fig. 5 of the main text, the ramp time was scaled 
with 14 as the tunneling rates and thus the band gap de¬ 
crease with larger I 4 . The switching between the lasers 
at 0.627r and 1.57r was done instantaneously within a 3 ms 
hold time between the ramps. It was confirmed experi¬ 
mentally that this switch does not lead to any measur¬ 
able excitation to other bands. This scheme was used for 
all measurements with the exception of the points with 
p < 1.57r in Fig. 2a of the main text, where a single laser 
was used with one ramp for p < ir and two ramps for 
p > TT, with a duration of 100 ms for each ramp. 


MODELS FOR FINITE PUMPING EFFICIENCY 

When transitions to other bands occur during the 
pumping process, the deflection of the cloud can be 
changed due to the different Chern numbers of the bands. 
For the measurements in the Wannier tunneling limit, 
we model the band occupations in a simple two-band 
model assuming that a small fraction of atoms is ex¬ 
cited non-adiabatically to the other band during each 
cycle. In the two-band model, these transitions will pre¬ 
dominantly take place when ramping over the symmet¬ 
ric double well configuration where the gap between the 
two lowest bands is smallest and atoms tunnel from one 
site to the next. This situation occurs twice within one 
pump cycle and we define the transfer efficiency as the 
fraction of atoms staying in the initially populated band 
during one ramp over a symmetric double well, which 
corresponds to half a pump cycle and thus one step in 
the COM position. 

The band populations after the m-th step can be ex¬ 
pressed as 


n(™)=e™.n(o) (S.ll) 

with 



Here, ni ( 712 ) denotes the fraction of atoms in the first 
(second) band, ei is the fraction of atoms excited from 
the first to the second band during one step and e 2 the 
one transferred from the second to the first. While the 
contribution from non-adiabatic transitions induced by 
the pumping is expected to be the same for ei and e 2 , 
the finite lifetime of atoms in the second band generally 
leads to e 2 > ei. 

When both ei and €2 are very small, one can expand 
e™ in powers of Ci and €2 and obtains to first order 

/1-mei me 2 \ (g^^2) 

mci 1 — me2 J 

Assuming a perfect preparation of the atoms in one band, 
i.e. = (1,0) or = (0,1), the resulting band 

occupations depend only on ei and € 2 , respectively. This 
model was used to fit both the measured COM position 
as well as the site-resolved band mapping data with a 
single free parameter. 

The expected displacement for both bands can be cal¬ 
culated by evaluating the COM evolution of the corre¬ 
sponding Wannier functions for one half of a pump cycle 
from (/? = 0 to TT. For each segment p € ](m — l)7r,TO7r], 
the combined motion is obtained by adding up the con¬ 
tributions from the two bands weighted with the respec¬ 
tive population fraction and After each half 

of a pump cycle (p = mir), the COM position can be 
expressed in terms of the Chern numbers of the bands: 

7 / m, \ 

^ iz • I sign(w)el^lw^°^ j (S.13) 

with u = (i^i, 1 / 2 )- 

For the site occupations, the even-odd distribution was 
determined by integrating the probability density of the 
Wannier functions over the even and odd site of the cor¬ 
responding double well. These were then weighted in the 
same way with and to obtain the fractions on 
even and odd sites as a function of ei and £ 2 . 

CENTER-OF-MASS POSITION IN THE 
PERPENDICULAR DIRECTION 

In addition to determining the COM position along 
the pumping direction, we also verified that the pumping 
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q) (2ir) 

Figure S3. Evolution of the COM position along the x- and 
j/-a:xis when pumping is performed along the i-axis. The grey 
points depict the COM displacement in the y-direction for the 
data points of Fig.2b in the main text, which are illustrated 
in red. The positions along y were evaluated using the same 
procedure as for the ones along x. Each point is averaged over 
ten data sets and the error bars depict the standard deviation. 


along X does not lead to any measurable displacement of 
the cloud along the perpendicular y-direction. For the 
measurement depicted in Fig. 2b of the main text, the 
COM positions along both axes are shown in Fig. S3. As 
expected, we do not observe any significant displacement 
of the cloud along y over multiple cycles in both pumping 
directions. 
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